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Transmission spectra of metallic films or membranes perforated by arrays of subwavelength slits 
or holes have been widely interpreted as resonance absorption by surface plasmon polaritons (SPPs). 
Alternative interpretations involving evanescent waves diffracted on the surface have also been pro- 
posed. These two approaches lead to divergent predictions for some surface wave properties. Using 
far-field interferometry, we have carried out a series of measurements on elementary one-dimensional 
(1-D) subwavelength structures with the aim of testing key properties of the surface waves and com- 
paring them to predictions of these two points of view. 



Early reports of transmission through arrays of sub- 
wavelength holes in thin films and membranes [3,0, 
enhanced well beyond conventional expectation |j] have 
motivated numerous attempts to explain the underly- 
ing physics of these surprising results. Since the early 
experiments were carried out on metal films, surface 
plasmon polariton resonances |^ |^ were invoked to ex- 
plain the anomalously high transmission and to sug- 
gest new types of photonic devices 0] . Other interpre- 
tations based on "dynamical diffraction" in periodic slit 
and hole arrays la| various kinds of resonant cav- 
ity modes in 1-D slits and slit arrays 0, 0| have also 
been proposed. Reassessment of the early claims by new 
numerical studies jlL] and new measurements [l^ have 
prompted a sharp downward revision of the enhanced 
transmission factor from ~ 1000 to ~ 10 and have moti- 
vated the development of an alternative approach based 
on a composite diffractive evanescent wave (CDEW) |0 . 
This model constructs a composite surface wave from the 
distribution of diffracted evanescent modes (the inhomo- 
geneous modes of the "angular spectrum representation" 
of wave fields 0|) originating at an abrupt surface dis- 
continuity such as a subwavelength-sized hole, slit, or 
groove when irradiated by a train of plane waves. The 
CDEW model exhibits three specific properties. First, 
the surface wave is considered to be a composite of modes 
labelled by the propagation vector component parallel to 
the surface, and evanescent in the direction normal to 
the surface. This composite "wave packet" exhibits well- 
defined, regular nodal positions spaced by a characteris- 
tic wavelength, Asurf ; second, the appearance of the first 
node at a distance of Asurf /2 from the structured edge 
implies an effective phase delay of 7r/2 with respect to 
the E-field of the external driving source; and third, an 
amplitude decreasing inversely with distance from the 



launch site with an overall effective range of a few mi- 
crons. 

We have fabricated 1-D structures (slits and grooves 
with subwavelength widths) in thin silver films deposited 
on fused silica substrates designed to test these features 
of the CDEW model. The optical response can be stud- 
ied with the structures facing toward (input-side experi- 
ments) or away from (output-side experiments) a distant 
coherent plane-wave source. Results from the input-side 
experiments, exhibiting light transmission interference as 
a function of slit-groove distance, have been reported 
elsewhere jl4j . We report here measurements of output- 
side, far-field intensity fringes arising from interference 
between propagating waves transmitted through the slit 
and surface waves launched at the slit but, after travelling 
along the surface, reconverted to outgoing waves at the 
groove (see Fig.EJ. Studies of fringe frequency, phase, 
and contrast as a function of slit-groove distance and 
groove depth provide new and complementary informa- 
tion to the previously reported input-side experiments. 

The subwavelength structures are fabricated by fo- 
cused ion beam (FIB) milling (FEI Nova-600 Dual-Beam 
system, Ga"*" ions, 30keV) in a 400 nm thick layer of sil- 
ver evaporated onto flat fused silica microscope slides. 
A low beam current (50 pA) is used in order to achieve 
surface features defined with a lateral precision on the 
order of 10 nm and characterised by near- vertical side- 
walls and a minimal amount of edge rounding. Since 
it enables delivery of a variable ion dose to each pixel 
of the writing field, FIB milling conveniently allows the 
multiple-depth topography characteristic of the present 
devices to be formed in a single, self-aligned step. A 2-D 
matrix of structures is milled into the silver layer. Each 
matrix consists of 63 structures, nine columns, separated 
by 1.5 mm, and seven rows, separated by 2 mm. The 
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FIG. 1: Goniometer setup for measuring far-field light inten- 
sity and angular distributions. See text for details. 

first column contains only slits with no flanking grooves. 
Light transmission through the slits in this column is used 
to normalise the transmission in the remaining columns. 
Variations in transmission through each of the elements 
in the "slits only" column provide a measure of the uni- 
formity of the FIB fabrication process. The remaining 
structures consist of slits flanked on one side by a groove. 
The groove-slit distance is systematically increased from 
the initial to final matrix positions. The square micro- 
scope slides themselves are 25 mm on a side and 1 mm 
thick. 

Measurements were carried out using a home-built go- 
niometer shown in Fig.^ Output from a diode laser 
source, temperature stabilised and frequency-locked to 
the Cs ^Si/2(F = 4) ^ ^P3/2(F = 4, 5) crossover feature 
in a saturated absorption cell, is modulated at 850 Hz by 
a mechanical chopper, injected into a monomode optical 
fibre, focused and linearly polarised (TM polarisation, H- 
field transverse to the slit) before impinging perpendicu- 
larly on the structure matrix mounted in a sample holder. 
The beam waist diameter and confocal parameter of the 
illuminating source are 300 /im and 33 cm, respectively. 
Throughout this series of measurements the laser power 
density was maintained ~ lWcm~^. The sample holder 
itself is fixed to a precision x-y translator, and individual 
slit-groove structures of the 2-D matrix are successively 
positioned at the laser beam waist. A photodiode detec- 
tor is mounted at the end of a 200 mm rigid arm that 
rotates about an axis passing through the centre of the 
sample holder. A stepper motor drives the arm at cali- 
brated angular increments of 2.05 mrad per step, and the 
overall angular resolution of the goniometer is ~ 4 mrad. 
The photodetector output current passes to a lock-in am- 
plifier referenced to the optical chopper wheel. Data are 
collected on a personal computer that also controls the 
goniometer drive. 

The structures consist of a single subwavelength slit. 



100 nm wide flanked by one subwavelength groove. Both 
groove and slit are 20 /xm long. We have carried out mea- 
surements with grooves of two different widths (100 nm 
and 415 nm) and depths varying from 32 nm to 256 nm. 
Here we report results only for the 100 nm grooves since 
the wider structures yield similar results. 

As illustrated in Fig.|21 the far-field intensity pattern 
should exhibit interference fringes between Et directly 
propagating through the slit and Eg radiating from the 
grooves after having been transmitted by the surface 
waves E'surf launched at the output side of the slit. The 
frequency and phase of the interference pattern is a func- 
tion of the optical path difference and any "intrinsic" 
phase shift (e.g. due to groove shape or surface wave 
phase lag) of Eg with respect to Et- Figure El shows in- 
terference fringes at three representative slit-groove dis- 
tances Xsg as a function of the goniometer detector angle 
9. The fractional surface wave amplitude a, normalised 
to the incoming plane wave amplitude Ej, a = Esmf/Ei, 
is estimated from the Kowarz model [lj| to be about 
95%. The remaining 5% constitutes the amplitude frac- 
tion S of the light directly transmitted through the slit 
Et- A further fraction (3 is reconverted to propagating 
light Eg at the groove site Xsg and interferes with Et- 
The intensity, Ig, of the superposed wavefronts can be 
expressed as 

Igie) = \Et + = \6E, + l3aE, ex^pijl^ (1) 

The normalised intensity Ig/Io can then be expressed, 
with r]o = a(3/5, 

^ cx 1 + + 2r/o cos 7 where 7 = /cqZo + (2) 
Jo 

with ^0 ~ sin 6* and ip — kgurfXsg + ^int- The re- 
lations between lQ,Xsg,0 are shown in Fig.|21 The fre- 
quency and phase of the interference pattern depend 
on the slit-groove distance through the terms ksnrfXsg, 
and kgXsg sm9. The term ipi^t represents the "intrinsic" 
phase shift due to groove geometry and surface wave gen- 
eration at the slit edge. With the goniometer detector 
oriented perpendicular to the structure plane {9 = 0), 
Eq.[21 simplifies to 

Y- OC l + 7]l+ 2?7o COs{ksuTfXsg + ifint) (3) 

Jo 

Figure 0] plots a series of measurements of the far-field 
intensity as a function of the slit-groove distance Xgg with 
the detector angle at 6* = 0. The interference term on the 
right side of Eq.O is fit to the data from which (p-mt can 
be determined by extrapolation of Xsg to zero distance. 
The magnitude of the surface wave propagation vector 
fcsurf = 27r/Asurf IS also determined from the fit. We 
mesure Asurt to be 811 ± 8nm in agreement with the 
input-side experiments reported earlier . The intrinsic 
phase ifint for this groove groove geometry (100 nm width 
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FIG. 2: Diagram showing interfering wavefronts, optical path 
difference between Et and -Eg, and far- field detection. 
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FIG. 3: Representative output-side interference fringes for 
structures with 100 nm wide grooves at three different dis- 
tances Xsg'. black, 0.543 /im; green, 1.845 ^m; blue, 4.991 /im. 
The phase value at zero observation angle (d = 0) is used to 
plot cos(f in Fig.2] 
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FIG. 4: Plot of the phase ksurfXsg + V'int as a function of 
Xsg with the detector oriented perpendicular to the struc- 
ture plane. Spectral analysis of the frequency spectrum of 
the fringes yields a determination of the surface wavelength 
Asurf = 811 ± 8nm. Extrapolation of the phase as Xsg ap- 
proaches zero, yields ipint ~ 0.32 n ± 0.02 n. 
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FIG. 5: Plot of rjo as a function of the slit-groove distance 
Xsg. The quantity rjo is fit by Eq.|S]with fitting parameters 
/i = 0.13 ± 0.01 and k = 0.12 ± O.Ol^m. 



and depth) is determined from the plot to be t^int = 

0.32 TTzb 0.02 TT. 

In addition to the frequency and phase of the interfer- 
ence we have studied the "visibility" or contrast of the 
output-side interference fringes as function of Xsg. The 
interference contrast is defined as 



C 



(4) 



where Imax, Imin are adjacent intensity maxima and min- 
ima of the fringes. According to Eq.|2]the contrast can 
be expressed as 
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Since ijo oc af3 — Eg/Ei, the fractional amplitude radiat- 
ing at a groove, a plot of r]o as a function of Xgg measures 



the dependence of this field amplitude (and therefore the 
surface wave amplitude) on the slit-groove distance. Fig- 
ure |S1 shows a plot of 77o as a function of Xgg for narrow- 
groove structures. The form of the fitted curve through 
the data points, an inverse distance dependence with an 
additive constant, is given by Eq.|Blwith fitting parame- 
ters /i, K as indicated in the captions of Fig.|Sl 



(6) 



rio{Xsg) = ^ ^J. 



^sg 



We can now compare these results to predictions of 
the CDEW model. Two key predictions are: (1) that the 
amplitude of the composite surface wave decreases as the 
inverse of the distance from the launch site, and (2) that 
there is an intrinsic phase delay (/?int of 7r/2 between 
the directly transmitted wave and i?surt the composite 



4 



evanescent wave. Figure HI shows that the contrast (and 
therefore the ampUtude of Eg) fits well a l/xsg behavior 
for about the first 3 microns, but then stabilises at a 
constant contrast. 

It appears therefore that there are two components to 
the surface wave amplitude: a rapidly decreasing com- 
ponent at short range followed by a constant component 
at longer range. Figure^shows that (fi^t extrapolates to 
0.32 7r± 0.02 (not 7r/2) as Xgg approaches zero. However, 
it is well known that grooves exhibit "organ-pipe" phase 
shifts and amplitude resonances when the effective depth 
is close to an integer number of quarter wavelengths [lOj . 
In order to investigate this contribution to the intrinsic 
phase we measured the contrast and phase as a function 
of groove depth. The results are shown in Fig. El The 
contrast indeed shows a maximum near 175 nm groove 
depth. Around this resonance the phase lag from the 
groove must be about modulo 2tt, and therefore any 
residual intrinsic phase between Et and Eg around the 
groove resonance must be attributed to the phase delay 
between the surface wave i?surf and the directly trans- 
mitted wave Ef . Figure [HI indeed shows that this resid- 
ual phase is close to 7r/2, the signature phase lag of the 
CDEW. Phase and amplitude data from the wide-groove 
studies jl6l | are consistent with results reported here. 

We conclude from these phase and amplitude results 
that the surface waves exhibit both CDEW-like and SPP- 
like properties. The initial decrease in interference con- 
trast, fitting well a 1/a; behavior, at small slit-groove 
distances is consistent with a diffractive surface pertur- 
bation at the slit edge. The persistance of essentially 
constant contrast at slit-groove distances greater than 
~ 2 fim indicates the presence of an an SPP-like long- 
lived mode propagating along the surface. It should be 
noted, however, that the expected wavelength of a pure 
SPP mode on a plane silver service is Aspp = 844 nm , 
but the measured (Fig.Q Asmf = 811 ±8. The reason for 
the discrepancy is unclear but this "output-side" deter- 
mination is consistent with the "input-side" results [l3| . 
When the "organ-pipe" phase contributions are taken 
into account, the results of Fig.lHI indicate that the the 
persistant SPP-like wave conserves the CDEW signature 
phase lag of 7r/2. The detailed nature of the conversion 
from a diffracted surface wave packet to an SPP guided 
wave has yet to be explained. 

Support from the Ministere delegue a I'Enseignement 
superieur et a la Recherche under the programme 
ACI-"Nanosciences-Nanotechnologies," the Region Midi- 
Pyrenees [SFC/CR 02/22], and FASTNet [HPRN-CT- 
2002-00304] EU Research Training Network, is grate- 
fully acknowledged as is support from the Caltech Kavli 
Nanoscience Institute and from the AFOSR under Plas- 
mon MURI FA9550-04-1-0434. Discussions and technical 
assistance from P. Lalanne, R. Mathevet, F. Kalkum, G. 
Derose, A. Scherer, D. Pacifici, J. Dionne, R. Walters and 



H. Atwater are also gratefully acknowledged. 




50 100 150 200 250 300 
Nominal groove depth (nm) 



FIG. 6; Intrinsic phase (fint (filled squares,left ordinate) and 
contrast (open circles, right ordinate) vs. groove depth. 
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